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Abstract 

A numerical value for the running electromagnetic coupling constant in the MS 
scheme is calculated at the low energy normalization scale equal to the r-lepton 
mass M T . This low energy boundary value is used for running the electro- 
magnetic coupling constant to larger scales where high precision experimental 
measurements can be performed. Particular scales of interest are the 6-quark 
mass for studying the T-resonance physics and the Z-boson mass Mz for high 
precision tests of the standard model and for the Higgs mass determination 
from radiative corrections. A numerical value for the running electromagnetic 
coupling constant at Mz in the on-shell renormalization scheme is also given. 



1 Introduction. 



Dimensional regularization (DR) and minimal subtraction (MS) are convenient and 
widely used technical tools for perturbative calculations in particle phenomenology 
H|, 0. In low orders of perturbation theory (PT) dimensional regularization does not 
give any decisive computational advantage. However, the high order PT many-loop 
calculations are rather involved and, in practice, only dimensional regularization sup- 
plemented by the recurrence relations based on the integration-by-part technique || 
allowed one to obtain new analytical results, e.g. [[§]. Minimal subtraction, being a 
simple method of renormalizing the dimensionally regularized PT diagrams, is now 
becoming a dominant way of theory parameterization in a form of the MS scheme Jj) . 
The renormalization in the MS scheme is mass independent that allows an efficient 
computation of renormalization group functions describing scaling of MS parame- 
ters. However, the mass independence of the renormalization procedure is physically 
inconvenient because decoupling of heavy particles is not automatic || . The phys- 
ical property of decoupling is restored within an effective theory approach with the 
explicit separation of different mass scales such that the parameters of neighboring 
effective theories (couplings, masses, ...) should be sewed (matched) near the point 
where a new scale appears. This machinery, worked out up to three-loop order in PT, 
allows one to compare theoretical results in the MS scheme for a variety of scales with 
a uniform control over the precision of PT calculation. In particular, this technique 
allows one to compare theoretical quantities extracted from the low energy data with 
results of the Z-boson peak analyses within the standard model (SM) of particle inter- 
actions. The high precision tests of SM at the Z-boson peak, completed at one-loop 
level, have shown a good agreement with theoretical results obtained from the low 
energy data. For the search of new physics and further tests of SM at the next level of 
precision the computations for many observables at the Z-boson peak should be done 



2 



with two-loop accuracy which presently is an actual calculational task. Because of 
the computational advantage of dimensional regularization in many-loop calculations, 
the high order PT results for theoretical amplitudes at the Z-boson peak tend to be 
obtained in terms of the MS scheme parameters which are natural quantities for the 
minimally-subtracted dimensionally-regularized diagrams. It was found that the use 
of the running electromagnetic (EM) coupling normalized at Mz in the MS scheme 
makes PT expansions near the Z-boson peak reliable and corrections small. However, 
contrary to the fine structure constant a, the running EM coupling in the MS scheme 
has no immediate physical meaning and its numerical value is not well known. At the 
same time, QED, being an old part of the SM, is well tested at low energies where 
the fine structure constant a is a natural interaction parameter defined in a physical 
manner by subtraction on the photon mass shell. The fine structure constant is very 
well known numerically that would make it a natural reference parameter for high 
precision tests of SM. However, because of a huge numerical difference between the 
values of the photon and Z-boson masses the use of the fine structure constant for 
calculations at the Z-boson peak generates large corrections in PT. For applications 
to high precision tests of the standard model with observables near the Z-boson peak 
0, one should transform a into a proper high energy parameter, i.e. into the electro- 
magnetic coupling constant at a scale of the Z-boson mass Mz (see e.g. || Then 
the large PT corrections are hidden (renormalized) into a numerical value of this new 
parameter which is more natural for describing the Z-boson peak observables than 
a. Therefore, a numerical value of the running EM coupling constant at Mz is a new 
important number which has been chosen as a standard reference parameter [K|. A 
difference of the numerical value for this parameter from a^ 1 = 137.036 . . . should 
be theoretically calculated. The change is accounted for through the renormalization 
group (RG) technique [11], [L^, [L3]. Because the fine structure constant is defined at 
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vanishing momentum it is an infrared sensitive quantity and a contribution of strong 
interactions into its RG evolution cannot be computed perturbatively: the infrared 
region is a domain of strong coupling that requires a nonperturbative (nonPT) treat- 
ment. The contribution of the infrared (IR) region is usually taken into account 
within a semi-phenomenological approximation through a dispersion relation with di- 
rect integration of low energy data. There has been a renewal of interest in a precise 
determination of the hadronic contribution into the electromagnetic coupling constant 
at Mz during the last years in connection with the constraints on the Higgs boson 
mass from radiative corrections in SM |14]. Some recent references giving a state-of- 



the-art analysis of this contribution are [jn| [U| [17], [18| . A quasi-analytical approach 
was used in ref. JL9|| where some references to earlier papers can be found (see also 



pHp). An extremely thorough data-based analysis is given in ref. |[22|| . However, 
the virtual lack of data for energies higher than 15 4- 20 GeV makes it unavoidable 
to use theoretical formulas in the dispersion relation. Fortunately, the theoretical 
results necessary for electromagnetic current correlators (the photon vacuum polar- 
ization function) are known in high orders of PT and are reliable at large energies. 
Therefore, the real value of dispersion relations is to find a boundary condition for the 
running EM coupling at a low energy normalization scale where data are accurate. 
If this low energy normalization scale is large enough for strong interaction PT to be 
applicable then the renormalization group can be used to run the initial value to any 
larger scale with very high precision. The running of the electromagnetic coupling 
constant can be defined in different ways depending on the renormalization proce- 
dure chosen. The evolution can be described in both on-shell and MS schemes: the 
corresponding /5-functions are available with high precision. The recent calculation 
of the numerical value for the running EM coupling at Mz with evolution in the MS 
scheme has been presented in ref. |S3 . 



4 



In the present paper I calculate a low energy boundary condition for the running 
electromagnetic coupling constant in the MS scheme using almost no experimental 
data but masses of ground states in the p- and yj-meson channels. The necessary 
IR modification of the light quark spectrum is determined by consistency with OPE. 
The theoretical parameters of the calculation are the strong coupling constant a s (M T ), 
the strange quark mass m s (M T ) and the gluon and quark vacuum condensates. The 
numerical values for these parameters accumulate a lot of information on low energy 
data presented in the standard rate R{s) of e + e~ annihilation into hadrons. Therefore, 
the present calculation compresses low energy data into the numerical values of several 
key theoretical parameters that allows one to perform an analysis of the IR domain 
necessary for determination of the low energy boundary value for the running EM 
coupling. The evolution to larger scales is straightforward and very precise within 
perturbation theory. 

2 Basic relations 

The relation between the running EM coupling constant in the MS scheme 

and the fine structure constant a is obtained by considering the photon vacuum 
polarization function. The correlator of the EM current 

J (Tj; m (x)jr(0))e^dx = (q,q u - g^q 2 )Tl # (q 2 ) (1) 

is defined with a generic scalar function U#(q 2 ). The particular scalar functions 
II(/i 2 ,g 2 ) and II os (g 2 ) are defined through the correlator of electromagnetic currents 
eq. ([I]) (and the generic function II # (g 2 )) but with different subtraction procedures to 
remove infinities. The first function U(p 2 ,q 2 ) is renormalized in the MS scheme and 
the second function U os (q 2 ) is renormalized by subtraction on the photon mass-shell 
q 2 = which implies the normalization condition I1 OS (0) = 0. Note that for the actual 
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calculation of II os (g 2 ) one can use dimensional regularization and the MS scheme in 
cases when IT(/i 2 , 0) exists, 



U os (q 2 ) = U(^,q 2 )-U(^,0). 
The relation between couplings and polarization functions in different schemes reads 

:]K + n( / Ag 2 ) = - + rw g 2 ). (2) 



a(n 2 ) ' a 

In the limit q 2 — > one finds 



3vr „. 9 . 37T 

+ n(/i 2 ,o) = — . (3) 



ct(/i 2 ) x ' a 

Eq. (H) is related to the Coulomb law for charged particles. For the potential of the 
EM interaction of two charged leptons one finds in the MS scheme 

_ 47m(/i 2 ) 1 

q 1 + -^nO 2 , q 2 ) 

This expression is \i independent because of RG invariance. Being expressed through 
the fine structure constant a the Coulomb potential reads 

2 4ira 1 

^ » = -IT i + jM) < 5 > 

with I1 OS (0) = 0. The limit of large distances 

4vra = - lim qV(q 2 ) (6) 

q 2 ^0 

gives the fine structure constant. In the Coulomb law eqs. (f|,|5]) q = (0, q) and 
q 2 = — q 2 . This makes q 2 in eq. (@) essentially Euclidean. We keep notation q 2 
for a positive number to stress the calculation in the Euclidean domain. Eq. (^) 
is just a relation between schemes in which the EM coupling is defined or finite 
renor malizat ion . 
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For massless quarks the limit q 2 — > in eq. (Q) requires care. The polariza- 
tion function n(// 2 ,0) cannot be calculated in PT if strong interactions are included 
because light quarks are essentially massless. 

Besides the MS running coupling constant a(/i), the on-shell running coupling 
a os (q 2 ), which can also be used in the Z-boson peak analyses, is defined through 

ex 

a os(<l 2 ) = : ~ ^ n , 2V «o S (0) = a. (7) 

1 T 3 7r J- L oslVl ) 

The numerical value for the on-shell running coupling a os (q 2 ) can be found from 
eq. (H) if a(/i) is known and IT(/i 2 , q 2 ) is calculable for a given q 2 . 

In the present paper I calculate the low energy boundary condition for the running 
EM coupling in the MS scheme, i.e. the value a(fio) at some fj,Q. A convenient scale 
is the r lepton mass M T which is large enough for strong interaction PT to work, 
i.e. Ho — M T . The value a(M T ) can then be run to other scales with the standard 
RG equation. The particular values of interest are m;, for T physics and Mz for high 
precision SM tests and Higgs boson search. The RG functions in the MS scheme are 
known with a very high accuracy that makes the running very precise numerically. 



3 Low energy normalization: formulas 

One needs a numerical value of the polarization function n(/i 2 , q 2 ) at q 2 = at some 
low normalization point /i 2 . There are lepton and quark contributions to the EM 
current (see a note about W bosons below). Because decoupling is not explicit we 
count only contributions of particles which are considered active for a given scale. 

3.1 Leptons 

For a lepton T with the pole mass Mi we retain masses that makes Il( / u 2 , 0) directly 
computable in low orders of PT where strong interactions are absent. The matching 
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condition reads 

nV,0)=b^-!-^-(--l^ln^-l 4 0(,r 



fi 2 a(/i 2 ) /45 3 /x 2 
M, 2 7T V 16 4 M? 



Note that (9 (a 2 ) corrections are also available []23| but they are totally negligible nu- 
merically for our purposes. With accuracy of order a there is no numerical difference 
between the fine structure constant a and the running coupling constant a(/z 2 ) in the 
RHS of eq. (§). For numerical estimates we substitute a. For fi = Mi the lepton '1' 
decouples completely in the leading order (which can be practical for the r lepton). 
Because the fine structure constant a is small numerically we do not resum the ex- 
pression in the RHS of eq. (§). Note that the expression ([8]) can basically be used at 
any fi. In a sense the matching for leptons can be done just at any scale of interest, 
for instance, at fi = Mz- For a lepton we use the pole mass Mi which is a meaningful 
parameter in finite order PT. Eq. (^ gives the leptonic part of finite renormalization 
between the running and fine structure constants in eq. (Q). In eq. (H) we neglect 
strong interactions (quark contributions) which appear in 0(a 2 ) order. If strong in- 
teractions are included then one cannot use PT with such a low scale as the electron 
or muon mass and the full IR analysis analogous to that done for light quarks (see 
below) is necessary. Eq. (Rp solves the lepton part of the normalization condition. 



3.2 Light quarks 

For the hadronic part of the vacuum polarization function we first consider a light 
quark contribution which is most complicated. For the light (massless) quarks the 
limit q 2 — ► in eq. (§) necessary to relate the running coupling to the fine structure 
constant cannot be performed in PT. This is, however, an infrared (IR) problem which 
is unsolved in QCD within PT. The low energy domain is not described in QCD with 
massless quarks by PT means and PT expressions should be modified for the limit 
q 2 — > in eq. (^]) to exist. Such a modification must not change an ultraviolet (UV) 
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structure of the correlators because RG invariance should be respected. Therefore, 
it is convenient to perform an IR modification using dispersion relations which give 
contributions of different energy ranges separately. There are three potentially IR 
dangerous quarks u, d, and s. For matching light quarks we work in rif = 3 effective 
theory, i.e. in QCD with three active light quarks. 

A note about notation is in order. We consider a generic light quark correlator 
normalized at the parton level to 1 (as for its asymptotic spectral density). Then 
we add necessary factors to account for color and/or charge structure. Thus, for u 
quark, for instance, 

n u (q 2 ) = N c e 2 u Ii li9h \^) (9) 

where e u = 2/3 is a w-quark EM charge and N c = 3 is a number of colors. For light 
quarks the PT part of the correlator is calculable for large q 2 and reads in the MS 
scheme (e.g. Ip5fl ) 

IT^ V, q 2 ) = 1 n ^ + I + a s (\ n ^ + § - 4 C (3)) 
9 /9 l2 u 2 /299 , x \ . u 2 34525 715 A/x 25,. A , , 

where a s = a s /ir, a s = a^(/i). Eq. (]T0|) is written for nj = 3 active light quarks. 
The limit q 2 —>■ cannot be performed because there is no scale for light quarks and 
no PT expression as eq. (|§|) is available. 

Still small momenta mean IR problems and we want to modify only the IR struc- 
ture of the correlator U h9ht (fi 2 , q 2 ). It is convenient to modify just the contribution of 
low energy states into the correlator which can be done through a dispersion relation. 
The dispersion relation reads 



p H 3 h \ s )ds 



IT^(q 2 ) = / F - (11) 







s + q 2 



where dimensional regularization is understood for p ll9ht (s). In fact, eq. ( |TTD can be 
used for bare quantities n z * 9ht (q 2 ) and p h9ht (s). The limit q 2 — > in the RHS of 
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eq. ( |Tl"D is IR-singular and cannot be performed if the PT expression for the spectral 



density p ll9ht {s) is used. Therefore, one should modify the low energy behavior of 
the spectrum where PT is not applicable. If such a modification is local (has only a 
finite support in energy variable s in eq. (|llD ) then it does not affect any UV prop- 
erties (fi 2 dependence) of U h9ht (n 2 , q 2 ) which are important for RG. The low energy 
modification is inspired by experiment: at low energies there is a well-pronounced 
bound state as a result of strong interaction. We, therefore, adopt a model of IR 
modification that the high energy tail of the integral in eq. (JTTJ) is computed in PT 
(duality arguments) that retains RG structure of the result while in the low energy 
domain there is a contribution of a single resonance. An IR modification is performed 
for contributions of u, d and s quarks. The massless u and d quarks interact with 
photons through isotopic combinations 1=1 (p-meson channel) and 1 = (cu-meson 
channel). For our purposes these two channels are completely degenerate and are 
treated simultaneously. The s-quark contribution is considered separately because of 
its nonvanishing (small) mass m s . 

For a generic light quark correlator in the massless PT approximation we introduce 
the IR modification 



where Fr, iur and Sq are IR parameters of the spectrum. Note that they are not nec- 
essarily immediate numbers from experiment. Substituting the IR modified spectrum 



P 



light 




light 




(12) 



( fL2|) into eq. ([□]) one finds 
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only input giving a scale to the problem. The IR modifying parameters Fr and so 
are fixed from duality arguments. 

Notice the difference in the 0(a 2 ) order between eq. fllPf) and eq. (|T3|) : in eq. (" p~3|) 
there is a new term — 37r 2 /8. This is so called V 2 ' correction (e.g. ||26|| ). It can be 
rewritten through ((2) = 7r 2 /6. 

To describe the IR structure of the correlator in representation given by eq. ( |13|) 
we use OPE with power corrections that semi-phenomenologically encode information 
about the low energy domain of the spectrum through the vacuum condensates of local 
gauge invariant operators [27]. The OPE for the light quark correlator reads 



n op V, q 2 ) = n'^ V, q 2 ) + ^ + o (^f \ . (w) 

The quantities (04,6) gi ye nonPT contributions of dimension-four and dimension-six 
vacuum condensates. These contributions are UV soft (they do not change short 
distance properties) and related to the IR modification of the spectrum. For the 
purposes of fixing the numerical values of the parameters Fr and Sq which describe 
the IR modification of the spectrum one needs only first two power corrections 1/q 2 
and 1/q 4 ; the coefficient of the 1/q 2 correction vanishes because there are no gauge 
invariant dimension- two operators in the massless limit. Computing the IR modified 
polarization function and comparing it with the OPE expression we find finite energy 



sum rules (FESR) for fixing the parameters Fr and so ||28|| . The system of sum rules 
has the form 

F R = so|l + a s + a 2 (/5oln^- + fc 1 +/5 )| + 0(a3), 
F R m 2 R = ||l + a s + a 2 ^ ln^ + A; 1 + ^|-(O 4 ) + O(a3). (15) 
Here (3 = 9/4 and 

2QQ 

h = - 9C(3) . 
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We treat {O4) as a small correction and take its coefficient function as a constant (the 
total contribution is RG invariant). Eqs. ([15]) fix F R and s through m R and (O4). 
Using higher order terms in OPE expansion ((Oq)/c[ 6 ) one can avoid substituting m R 



from experiment because within the IR modification given in eq. (|12|) the IR scale is 



determined by the dimension-six vacuum condensate {Oq) [28|. We do not do that 
because the primary purpose is to find the normalization for the EM coupling and not 
to describe the spectrum in the low energy domain. The use of the experimental value 
for the resonance mass m R makes the calculation more precise because the numerical 
value for (0 6 ) condensate is not known well (cf ref. [f29| ). 

The leading order solution to eqs. (|l|) (upon neglecting the PT and nonPT cor- 
rections) is given by the partonic model result Sq = 2m R , F R = s = 2m R , which 
is rather precise. This solution has been used for predicting masses and residues of 
the radial excitations of vector mesons within the local duality approach when the 
experimental spectrum is approximated by a sequence of infinitely narrow resonances 
|30| . Such an approximation for the experimental spectrum is justified by theoretical 



considerations in the large N c limit |JT|] and by the exact solution for two dimensional 
QCD ||32|| . For the experimental spectrum of infinitely narrow resonances the local 
duality approach means averaging over the energy interval around a single resonance 
pOfl . It is expected to be less precise than the global duality method in which the 
average is performed over the entire spectrum. However, within the global duality 
approach only the total contribution of all hadronic states can be studied while the 
local duality can be used even for first resonances and can predict characteristics of 
individual hadronic states. 

An accurate treatment of eqs. ( |15|) gives the solution 

.2 



*l = 2 l {1 + a , + al ^jL + h + l^ + ^. (16) 
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In the solution given in eq. ([16]) only linear terms in the nonPT correction (O4) are 
retained. This is well justified numerically. The a^O/i) terms are neglected because 
the coefficient function of the (O4) condensate is not known with such precision. In 
eq. (|l3| ) the scale parameter is Sq while we solve the system (|l(j) in terms of m# 
that we identify with the resonance mass and take from experiment. Therefore, we 
express the PT scale sq through itlr according to the solution given in eqs. (|16D . The 
expansion of the log-term in eq. fll3|) reads 



/i 2 u 2 ft 2 (0 4 > 

— = m 7, a e -H 

s 2m\ 2 s 2m% 



With these results one finds an expression for the IR modified polarization function 
of light quarks at the origin 

nX orf (/A 0) = 2 1 1 + a s + a 2 (ft In ^ + k x + -ft) } + ^ 

f A) 2 <04> h .,5 / /i 2 (CM 55 \ 

2 /9 l2 /i 2 /299 , ,\ /i 2 34525 715 _ 25 A . r . 3tt 2 \ 

+a 2 - In 2 + 9C(3) In -*-=- + C(3) + — C(5) • 

s V 8 2m l ^24 sv 7 2m\ 864 18 SV ; 3 SV ; 8 / 

Here the first line gives the resonance contribution while the rest is the high energy 
tail (continuum contribution) which is computed in PT. Finally, 

HSU/A 0) = 2 1 1 + a s + a 2 (ft In £^ + k x + 2ft) } + g| 

+ ln JfL_^ a 2 + 5 f ln j4 + --4C(3)^l 
2m| 2 s 3 S V 2m\ 12 sv 7 

,2 



ft, 2 /i 2 , , /i 2 34525 715,, , 25 . . 3tt' 



+ a 2 I — In" + h In + - — C(3) + — C(5) - — ] . (17) 

s V 2 2m 2 R 2m 2 R 864 18 S V ; 3 S ' 



Eq. ( |TTD gives n / i ^ mod (/i 2 , 0) as an explicit function of the nonPT scale m# (to be 
taken from experiment) and theoretical quantities a s and (O4). The choice of the 
numerical value for a s is discussed in detail later. 
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The condensate of dimension-four operators for light quarks is given by 

(0 4 > = y( 1 + l as ) { V° 2) + 21x2 i 1 + l as ) (m " + md) ^ Uu) + ^ d)) ■ (18) 

We retain small corrections proportional to the light quark masses and treat them in 
the approximation of isotopic symmetry for the light quark condensates (uu) = (dd) 
which is rather precise for u and d quarks. The quark condensate part of eq. ( plf ) is 
given by the PCAC relation for the ir meson 

(m u + m d )(uu + dd) = -f%ml . 

Here f n = 133 MeV is a charged pion decay constant and m n = 139.6 MeV is 
a charged pion mass. For the standard numerical value of the gluon condensate 
(^G 2 ) = 0.012 GeV 4 [0 and a s = 0.1 one finds 

71-2 ^ 7 \ ,Ois r, 2 /i , 1 „ \ {2 2 _ n n97 n^T± 



(o,) = y{ 1 + e a V V } ~ V + 3 a V ^ ' = °' 037 • (19) 

For the most important contribution of u and d quarks (the u-quark contribution 
is enhanced by factor 4 because of its doubled electric charge in comparison to the 
other light quarks) the relation sq = 2m 2 p , where m p = 768.5 MeV is a mass of 
the lowest (p meson) resonance in the non-strange isotopic 1=1 vector channel, is 
rather precise numerically. The gluon condensate gives a small correction to the basic 
duality relation for light quarks Sq = 2m 2 R . Note that we do not identify Fr with 
the experimental number available from the analysis of the p-meson leptonic width 
but treat it as an IR modifying parameter which should be fixed from the consistency 
requirement with OPE. It is close numerically to its experimental counterpart because 
it is known that OPE gives rather an accurate description of the physical spectrum 
if vacuum condensates are included. In the present paper we stick to a theoretical 
description of the IR domain and use the lowest resonance mass as the only input 
for the IR modification. The same is true for the 1 = channel where the lowest 
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resonance is the uj meson with a mass m w = 781.94 MeV. We do not distinguish 
these two channels. We consider parameters F R and So as the IR modifiers fixed 
theoretically through OPE and do not attempt to substitute them from experiment 
(using leptonic decay widths for F R or the shape of the spectrum for s ). 

Note that the IR parameters of the spectrum F R , m R and s are \i independent. 
It can be seen explicitly from eqs. (|15|) . 

The n/ = 3 effective theory is valid only up to q 2 ~ m 2 and, formally, there are 
corrections of order q 2 /m 2 []31|. However, in the case of current correlators these 
corrections are small [34, 35]. 

For the s quark there are also corrections due to m s which change slightly the shape 
of the spectrum and the consistency equations for the IR modifiers. We consider m s 
as an additional IR modifier which does not affect UV properties (renormalization in 
the MS scheme is mass independent) and treat it as a power correction. We write 
OPE for the s quark in the form 

n OPB ' V, q 2 ) = n^ 2 , q 2 ) - ^$ + M + o ( ^f] . 

The system of equations for fixing the parameters F Rs and Sq s reads 



F Rs + 6m 2 s = s 0s {l + a s + a 2 (/3 ln ^- + h + (3 )\ + 0(a 3 s ) , 

s 0s ) 



Here 



F s m 2 Rs = S fU + a s + a 2 s U\n^ + k 1 + ^)\-(Ot) + 0(a 3 s ). (20) 



(Ol) = y (l + \a) (yG 2 ) + Stt 2 (l + ±a s ) m s (ss) 



is a dimension-four contribution in the strange channel. One finds the solution to 
eqs. (EC]) in the form 



s , = 2m R s 1 + - a s +— —(l-a s )-6m 

\ £ I lib 



Rs 

^ = 2{l + a s + a 2 Lln-^ + A ;i + ^) + ^-12^i. (21) 
m 2 Rs { V 2m Rs 2 /J m Rs m 2 Rs 
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The correction due to m 2 is not large. Instead of eq. ([T7|) one has 



nSKkV, o) = nSU/A o) - (22) 



2 



m 



Rs 



and rriR S = and should be used in the first term of eq. (|22) instead of m p and 
(O4). Here m,, = 1019.4 MeV is a mass of the yj-meson which is the lowest resonance 
in the strange channel. A numerical value for (Ol) is obtained as follows. We use the 



relation (e.g. f36|1 ) 

2m, 



m u + m d 



25.0 



and the phenomenological result (ss) = (0.8 ± 0.2) (mm) |37j to find 

m s (ss) = 12.5 ■ 0.8 ■ (m u + m d )(uu) = -5.0 ■ f 2 ml = -0.0017 GeV 4 . (23) 

One could also use the PCAC relation for the K meson 

(m s + m u )(ss + uu) = -fKm 2 K + 0(m 2 s ) 

with fx = 1-17/tt and = 493.7 MeV. Note that the PCAC relation in the strange 
channel is valid only up to terms of order m\ which are not completely negligible 
numerically compared to the pion case [SB]. Therefore, we use the result given in 



eq. (|2§). For the standard value (^G 2 ) = 0.012 GeV 4 [0 and a s = 0.1 one finds 

(Ol) = y (l + 7 -a s ^j (^G 2 ) + Stt 2 (l + ia s ) (-5.0/> 2 ) = -0.0965 GeV 4 . (24) 

The correction due to m s (ss) is dominant in the dimension four contribution in the 
strange case. Because the dimension-four terms represent only small corrections to 
the leading results for the correlators in eqs. ( jr7|J22] ), the precision with which they 
are calculated suffices for our purposes. 

For the absolute value of m s to be substituted into m 2 correction we use the results 
of recent analyses [^] and take m s (M T ) = 130 ± 27 exp ± 9 th MeV. For m Rs = = 



1019.4 MeV one finds 

,2 



^ = 0.0163 
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which is a small expansion parameter that justifies the treatment of m 2 s contribution 
as a correction. 

Note that there are attempts to use constituent masses for the light quarks and to 
estimate the polarization functions in the way it is done for leptons or heavy quarks. 
Besides being ad hoc (and not supported by experiment) this IR modification of 
the light quark correlators contradicts OPE and/or the local duality over the energy 
interval from the origin to 1 4- 2 GeV. 

Thus, eqs. (|T^,|T^) represent a semi-phenomenological subtraction for a light quark 
correlator at q 2 = based on the IR modification of the spectrum consistent with 
OPE. Some mismatch with OPE in orders higher than 0(l/q 4 ) which is possible 
because of simplicity of the IR modification is neglected. It is justified because we 
need only integral characteristics of the spectrum for calculating n/^ od (/i 2 , 0) and 
are not interested in the point- wise behavior of the spectral function P/j^o^s) which 
is used as an auxiliary quantity in this particular instance. 

3.3 Heavy quarks 

Matching heavy quarks is straightforward and is similar to that of leptons. It is done 
within PT. For a heavy quark 'q' with the pole mass m q ^> Aqcd one has 

rF(/i 2 ,0) = N c e 2 q U heavy {fi 2 ,0) 

where U heavy (fi 2 , 0) is a generic contribution of a heavy quark to the polarization 
function KI 



n-V, 0) = In 4 + e£ (« + \ In + a, (H + In 

m l H 7r y 16 4 m z q J \ 4 m l ) 

2 /41219 917 / 4. n 2 \ 607 Wo , 
+ °- (ISM " 1296"' + ( 4 + 3 1,1 2 - 3"') C(2) + 144 C(3) 

+ (I* - ?) .n 4 + f H - l, tl ) to' 41 1 + O ia W.) ■ (25) 



36 9 "7 "' m\ V24 12 7 m\ t 
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Here ni is the number of quarks that are lighter than a heavy one, a s = a s /ir is the 
strong coupling constant in the effective theory with n\ + 1 active quarks normalized 



at the scale \x. Numbers in eq. (|25|) are given for the pole mass of a heavy quark. 
We neglect the (known) EM contribution of order a 2 because it is smaller than the 
unknown term of order a^. Eq. (p5|) gives a contribution of the corrected partonic 
model, i.e. that with a heavy quark loop in the leading approximation. There is also 
a contribution of heavy quark loops to the light quark vacuum polarization function 
that should be taken into account in constructing the effective theory with a decoupled 
heavy quark. This contribution is small. It reads El 



n— o) = (± 4) {—T^k-T^hj- (26) 



Eqs. ( p5| , p6|) are used for c and b quarks. Note that these formulas cannot be used 
for s quark. Indeed, because of a s corrections the PT scale in eq. (^) is effectively 
equal to m q and is too small for PT to be applicable in the case of the strange quark 
since m s ~ Aq C d- 



4 Low energy normalization: numerics 

In previous sections the necessary contributions due to fermions have been written 
down. We are not going to consider scales larger than Mz therefore bosonic contri- 
butions into the EM current and polarization function (namely, W boson loops) are 
not taken into account. The above equations describe an effective theory without W 
bosons which decouple at energies smaller than M z and should be added separately 
for the Z-boson peak tests. 

A numerical value of the strong coupling at low energies is rather important for the 
whole analysis. The estimates of the strong coupling numerical value at low scales are 
usually based on the r lepton decay data. Within a contour resummation technique 
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42|, |H the value obtained is a^(M?) = 0.343 ± 0.009 ex . p . Within a RG invariant 
treatment of ref. @ a slightly different value otf\M*) = 0.318 ± 0.006 exp ± 0.016 t/l 
has recently been found. The uncertainty is due to the experimental error and due 
to truncation of the series which is estimated within an optimistic scenario that 
higher order terms are still perturbative (no explicit asymptotic growth). Note that 
even for the optimistic scenario with a reduced theoretical error as compared to the 
conservative estimates, the theoretical error dominates the total uncertainty of the 
coupling. Contour improved results include a special resummation procedure for 
treating contributions generated by the running which does not necessarily improve 
results but definitely changes them in comparison with the finite order estimates at 
the present level of precision. The change is still within the error bars that makes two 
procedures consistent. We use the value af\M^) = 0.318 ± 0.017 as our basic input 
for the low energy strong coupling. The central value a^(M^) = 0.318 corresponds to 
af\Mz) = 0.118 when it is run with a four- loop /3-function and three- loop matching 
at m c and rut, thresholds. 

Now one has everything for numerical analysis. First the value of the running EM 
coupling constant computed in rif = 4 effective theory at [i = M T that is a convenient 
normalization point at the physical mass scale is given. Note that the c-quark pole 
mass is rather close to this value. In fact, the recent estimate is m c = 1.8 ± 0.2 GeV 
and we take m c = M T = 1.777 GeV as a central value, i.e. m c = M T ±0.2 GeV. Thus, 
the low energy normalization value a^(M T ) is computed with three active leptons 
and four active quarks. 

For the lepton contribution we use the lepton masses M e = 0.5110 MeV, = 
105.66 MeV, M T = 1777 MeV 0. These values are extremely precise therefore we 
use them as exact and assign no errors to them. We neglect the difference between 
the running EM coupling a and fine structure constant a in corrections (which results 
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in 0(a 2 ) shift that is numerically negligible). We use a 1 = 137.036. According to 
eq. (H|) leptons give 

A^(M 2 ) = £ rf(M T 2 ,0) = 1 + ln-£ + ln-| + ln-£ + 



47ry v My m? y i6 7T 

= 21.953 + 0.058 = 22.011 (27) 

where the first number is obtained in the limit a = 0. The a correction is almost neg- 
ligible for the normalization at the scale M T . Note that r lepton gives no logarithmic 
contribution at the scale /i = M T . 

The 0(a 2 ) correction for the lepton contribution in the MS scheme is also available 



24| . This correction is parametrically small and there are no unexpectedly large 
numerical coefficients (in fact, they are also small) that makes the parametric estimate 
based on the counting of powers of a rather precise. The sum of contributions of three 
leptons in 0(a 2 ) order is completely negligible and we treat the leptonic contribution 
in eqs. (f|^7]) as exact. 

Light quarks. From eq. (|i~7D with m R = m p , and rriR S = m v one finds for the total 



light quark contribution A (M 2 ) 

A uds (M 2 ) = A U (M 2 ) + A d (M 2 ) + A S (M 2 ) = A P {M 2 ) + A"(M 2 ) + A^(M 2 ) 

_ ^A light {M 2 ) + -A li9ht (M 2 ) + -A li9ht ~ s (M 2 ) = -A li9ht (M 2 ) + -A li9ht ~ s (M 2 ) 
3 3 3 3 3 

= 9.13662 + 5.32853a s + 24.9086a 2 

= 9.13662 + 0.53937 + 0.255214 = 9.9312 . 

Because the calculation is explicit we can give the previous result in more detail 
showing all different contributions 

A uds (M 2 ) = 9.11165 + 0.539367 (-^—} + 0.2552 ( 

v T> vo.ioiy vo.ioiy 

+ 0.08865 ( ^—s) - 0.0488 ( ^ - V + 0.0149 ( ^ j) . (28) 

V0.037 GeV 4 y V130 MeV J V - 0965 GeV / 
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The IR part of the spectrum (resonances) and the partonic quark approximation give 
a dominant contribution. The QCD perturbative corrections and power corrections 
due to m s and (Of) condensates are small. The error is 

6A uds (M 2 T ) = 10.55a, + 0.09^ - 0.1^ - 0.015^ . (29) 

(04/ m s (CI) 

Variations 5 {Of) and 5{0^) are not completely independent - both quantities contain 
a variation of the gluon condensate. Also the error of a s and that of the gluon 
condensate are correlated (see, for instance, |5J). For estimating the total error of 
A uds (M 2 ) through less correlated quantities one could rewrite power corrections in 
eq. ( P%D in the basis of the gluon and strange quark condensates [[K]]. Because the 
correlation is not well established numerically we neglect this effect. We consider the 
errors of the strong coupling a s , of the gluon condensate for 6(04), of the strange 
quark mass m s , and of the strange quark condensate (ss) for S(Of) as independent 
and use 5a s = 0.017/% = 0.0054, 8(0^) / (O4) = 1/2 due to the gluon condensate, 
5m s /m s = 0.28, 5(C|)/((9|) = 1/4 due to (ss). With these (conservative) estimates 
of uncertainties one finds 

5A uds (M 2 T ) = ±0.057| as ± 0.045| m ± 0.028| ms ± 0.004| {O|) . 

The dominant error is due to 5a s . The gluon condensate gives a sizable error because 
it is enhanced by the charge structure of light (mainly u) quarks and because its 
uncertainty is taken to be very conservative to compensate for the possible correlation 
with a s . The strange channel is suppressed by factor 1/3 in the total sum of light 
quark contributions and its specific features only slightly affect the result: in the rest 
it is quite degenerate with u and d channels. The total error for the light quark 
contributions added in quadrature reads 

5A uds (M*) = ±0.078. 
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The final result for the contribution of light quarks into the low energy normalization 
of the running EM coupling is 

A uds (M 2 ) = 9.9312 ± 0.078 . (30) 

We retain some additional digits at intermediate stages just for computational pur- 
poses. 



For the c quark we use eqs. (|25|j26|) . The strong coupling constant in rif = 4 
effective theory is found by matching between strong coupling constant in n/ = 3 and 
n/ = 4 effective theories. 



Matching at the pole mass scale mp for the strong coupling has the form 

a^\m 2 P ) = a^ +1 \m 2 P ) (l + C 2 a^ +l \m P f + C 3 a^ +1 \m 2 P ) 3 + 0(a 4 s )) (31) 
where 

C 2 = ~^, (32) 

80507 A/ , 2 A/ W1 , 58933 ni A, , 2479\ , , 

C 3 = C 3 - -( h +3 + — C 2 + . 33 

3 27648 SV ; 9 SV A ; 124416 9 V V ; 3456/ V ; 

We solve (inverse) eq. (^Tj) perturbatively and find the expression 

a^ +1 \m 2 P ) = 4 ni \m 2 P ) {l - C 2 a { s ni \m 2 P ) 2 - C 3 a^\m 2 P ) 3 } (34) 

which is used for determination of the couplings in neighboring effective theories at 
their boundary scale that is chosen to be the pole mass of the heavy quark. Matching 
at m c = M T = 1.777 GeV (we remind the reader that the numerical value of the 
c-quark mass is chosen to be m c = M T ± 0.2 GeV) with af'(M 2 ) = 0.318 gives 
a^\m 2 c = M 2 ) = 0.102 or a^\m 2 c = M 2 ) = 0.320. This value for the strong 
coupling is used in eq. ( |25|) for the calculation of the c-quark contribution to the finite 
renormalization of the EM coupling. Note that though one computes with af\M 2 ) 
it can well be identified numerically with af\M 2 ): the change due to matching is 
tiny and is much smaller than the error of af\M 2 ). 
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According to eqs. (|25| , p6|) we have 

A c (M T 2 ) = n c (/i 2 = M r 2 ,0) 
= 0.00387 + 0.00474 + 0.51001 + 0.32817 = 0.84679 (35) 



where the first term is EM contribution, the second one is loop contribution (eq. 
and the last two terms give PT expansion of direct contribution (eq. (|25|)). One sees 
that the EM and loop contributions are much smaller than the direct contribution. 
Convergence of PT series for the direct contribution is not fast though. 

The uncertainty of the c-quark contribution is straightforward to estimate in view 
of explicit formulas. The main error comes from the uncertainty of the c-quark mass. 
In the next-to-leading order one has from eq. (|25f) 



5A C (M 2 T ) = 4(1 + a.)— = + *.)— = ±0-330 (36) 

3 m c 3 m c 

for m c = M T ± 0.2 GeV and a s = 0.1. This is a very large uncertainty. The con- 
tribution A C (M 2 ) in eq. ( |3"5] ) is small because the c quark almost decouples but the 
uncertainty of A C (M 2 ) is large. The uncertainty is mainly given by the variation of 
ln(M 2 /m 2 ) in eq. (KH) and is independent of an absolute value of the contribution. 
For the central value m c = M T one would find a vanishing contribution in the leading 
order but its uncertainty would almost stay unchanged and equal to 0.330. Also the 
c-quark mass is not very large and the convergence of the PT expansion in eqs. (^5],[35D 
is slow. 

Note that for estimating the c-quark contribution we do not take into account 
the a s uncertainty. The reason is that uncertainties of the quark mass and of a s are 
strongly correlated. Indeed, to the leading order one can find from eq. (|25|) 

Sn heav y(M*,0) = —8a s (37) 
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for independent variation of a s . However, eq. fl25|) can be rewritten in terms of the 
running mass rh c (fi 2 ). To the first order in a s the relation between masses reads 

m c = rh c (^)\l + a s (fi 2 ) ( In — + f]\ (38) 



m 2 (jj 2 ) 3 



that leads to the change in eq. 



In At + a s | — + In At I -> In _ n r ON + a s \ — - In ■ 



m 2 \4 m 2 / m 2 (/i. 2 ) y 12 m 2 (/i 2 )y 

The NLO result for the polarization function in terms of the running mass now reads 

n^ s (^,0) = ln-^ + a s f^-ln. ''" ^ 



m c(/ i2 ) \12 m 2 (/i 2 ) J 
This result leads to the uncertainty 

/ 1 3 M 2 \ 

^ h rZ2 ass (M 2 T , 0) = - In J K 

which is by factor 4 smaller numerically than the previous result eq. (|3~?D. The rest 
of the uncertainty is now in the relation between the pole and running mass given 
in eq. (|38|) that represents a regular change of variables in the finite order PT and 
is under rather a strict control. We work with the pole mass and assume that the 
uncertainty of the polarization function at the origin is saturated by the uncertainty 
of the pole mass. It is also assumed that the uncertainty of the pole mass is estimated 
such that it includes an uncertainty of a s . 

The total finite renormalization between the fine structure constant and the MS- 
scheme EM coupling at M T is given by 

A (4) (M 2 ) = A lept (M 2 ) + A uds (M 2 ) + A C (M 2 ) 

= 22.0109 + 9.9312 + 0.84679 = 32.7889 (39) 



that leads to 



3vr 37r (A \ „. 37r . 

A (4) (M 2 ) = 32.7889 . (40) 



«( 4 )(M 2 ) a y T ' a 
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The low energy normalization value for the EM coupling in the MS scheme reads 

1 



133.557. (41) 



a( 4 )(M T 2 ) 

We now consider the uncertainty of this central result. The lepton contributions are 



treated as exact so the number from eq. (27]) has no errors. The errors due to light 
quarks are given in eq. (0). Note that one could reduce the sensitivity of A uds (M 2 ) 
to a s , the error of which dominates the total error in eq. (0), by taking Fr from 
experiment through the leptonic decay width of the p meson (and of the u and tp 
mesons in other light quark channels). Then the resonance contribution Fn/m R does 
not depend on a s . The first duality relation fixes Sq immediately using the fact 
that power corrections of order 1/q 2 are absent in the OPE. However, this procedure 
introduces an experimental error due to an uncertainty in the p-meson leptonic decay 
width 

T p ee = 6.77 ±0.32 keV. 

This uncertainty leads to almost the same error for the final quantity A uds (M 2 ) as 
the uncertainty in a s . It seems to be natural. Indeed, the strong coupling at low 
energies is extracted from the r data in which the p-meson contribution constitutes 
an essential part. This example shows how the coupling constant encodes information 
on the experimental data. Another point about using T p ee for the lowest resonance 
contribution is that the consistency with OPE is less strict for such a procedure 
(no dimension-four operators participate). Still having in mind the possibility of 
further improvement through experiment we consider our estimate of the error given 
in eq. (^) as rather conservative. 

The uncertainty of the c-quark contribution is given in eq. (^§). Collecting all 
together one finds the final prediction 

A {4) (M T 2 ) = 32.7889 ± 0.078^ ± 0.330 c (42) 
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and 



l7r ¥L _ A W(M T 2 ) = — - (32.7889 ± 0.078 %w ± 0.330 c ) . (43) 



a( 4 )(M r 2 ) a v Tl a 

Eq. ( ^|) is the main result for the low energy normalization of the running EM 
coupling. For the coupling itself it reads 

^(4)^2) = 133 ' 557 ± °- 0083 ^ ± 0.0350 c (44) 

and 

«( 4 )(M r 2 ) = 1.0261a. 

This value a^(M^) (or equivalently A^(M 2 )) represents the boundary (initial) con- 
dition for the running. With this value known the EM coupling can be run to other 
scales. The final goal is Mz = 91.187 GeV where high precision tests of SM are done. 
As will be seen later the running itself is very precise numerically and the main uncer- 
tainty in the running EM coupling at larger scales is due to the boundary condition 



eq. (f44|). The boundary condition eq. ( f44"|) has rather a big uncertainty because of 
the error of the c-quark mass mainly. The uncertainty due to the light quark contri- 
bution is reasonably small. It is dominated by the error in a s (M T ) which is mainly 
theoretical, i.e. related to the truncation of PT series used for describing the r-lepton 
decay data. The uncertainty in a s (M T ) can be reduced if some other sources for its 
determination are used in addition to the r system. Reducing the c-quark pole mass 
uncertainty requires more accurate treatment of the threshold region of cc production 
which is rather a challenging problem in QCD. 

5 RG evolution in MS scheme 

With the boundary value known at sufficiently large scale one can perturbatively run 
the EM coupling to larger scales. The final goal is the determination of the numerical 
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value for the EM coupling at Mz where high precision tests of the standard model 
are performed. The running itself (as a functional) is extremely precise because (3- 
functions are very well known. The precision of running is affected by the initial value 
of a s which is chosen to be af\M^) and by the 6-quark mass m^. We discuss them 
in detail later. 



5.1 Basic relations for RG evolution 

For the evolution between the r-lepton mass M T = 1.777 GeV (numerically m c = M T ) 
and Mz = 91.187 GeV the number of active quarks is either 4 or 5 and only one 
threshold at mj is encountered. The evolution equation (running) is written in the 
form 

2 d f 3n \ / 3a\ {10 l n \ (\1 l.\a 

"*? (spy ) = 3 1 1 + **) + (t + 5*) + (is + s 9 ') n 

- (| + ±h) + a.ft Q0D (a,) • (45) 
Here Q\, is a parameter for the 6-quark presence, rif = 4 + 9b- From M T to one has 



rif — 4 and 9b = while from m& to Mz one has rif = 5 and 0& = 1. In eq. (45) the 
strong coupling a s (/i 2 ) obeys RGE 

^a s (fi 2 ) = ^a s (fi 2 )) + al^(j:e^ (46) 

with 

/3(a.) = -a 2 s ((3 + faa. + foa 2 s + (3 3 a 3 s ) + 0(a 6 s ) (47) 
being the strong interaction /^-function. In QCD one has 



3 3 7 1 V144 72 

,2/38551 7595 77 2 55 \1 
1 15552 ^ 7776 06 ~ 3888^ " 54 C(3)(1 + ^ b) J I 



_, /2 l„\ 2 /55 5 
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where the first two lines give the 'direct' contribution and the third line gives light-by- 
light contribution which is written separately because of its different color structure. 
This result is obtained from the photon renormalization constant given in H and 



explicitly written in jy], [18). It was used in ref. |23| for calculation of the evolution 



of the EM coupling constant. Numerically, one finds 

h QCB (a s ) = (y + 3^) (l + a.(1.993 - O.1530 6 ) + a 2 8 (1.26 - 2.209 b - 0.020^ 

+ ^(-0.55 + 0.55^-0.14^). (49) 



Coefficients of the EM /3-function in eqs. (^,|^,^) are very small that makes the 
convergence of PT series for the evolution very fast. Eqs. (45,46) should be solved 
simultaneously. However, the EM coupling a(/i) is small, therefore, we neglect its 
running in corrections and substitute there the value numerically equal to the fine 
structure constant a. Then one has to integrate the trajectory of the strong coupling 
a s (n) which is given by the solution to RGE (P?p. The a correction in the strong 
coupling /5-function is numerically of order a 2 s and formally should be retained if 
terms in the /3-function are retained. However, the main contribution to the running 
is given by the partonic part of the EM /3-function in eq. (flSP, i.e. independent of 
both EM and strong couplings. Other terms give only small corrections. As for 
practical calculations, one can do everything numerically, however, it happens that 
the two-loop running gives almost the same result as the exact treatment. With the 
two-loop accuracy the integration can be done analytically in a simple form. Indeed, 
for P(a s ) = —f3oa 2 s — (3\cP s one finds 

/ a s {£)d In f = —In , 

U a ^J = -J^K^^M)) (50) 
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where the NLO solution for the running coupling a s (/z) is given by 

In^j =*(«.) = / _ m + ft{) = ^ + /f ln l^ftj- (51) 
In NNLO it is also possible to perform integration explicitly but results are too 



awkward to present here. In fact, the NLO integration as given in eqs. (|50| , |5ll) is rather 
precise numerically and can be used for preliminary estimates. We, however, avoid 
any approximation of this sort (cf ref. |23|) and give numbers for a direct numerical 
treatment of RG equations (|5|,[|(|) with the four-loop strong coupling /^-function from 
eq. (47) and /i^ CD -function from eq. (|48|). 



The solution to RGE can be used for the range of fi where the corresponding 
effective theory (with a given number of active leptons and quarks) is valid. Because 
decoupling is not automatic one should explicitly take into account thresholds. 

5.2 Running to mj 

The first scale of interest is m& where there is an important physics due to bb produc- 
tion near threshold and accurate data on T resonances (note that the real threshold 
energy is, in fact, 2rrib but the matching is defined at m&). We use m& = 4.8±0.2 GeV 
as determined in ref. p9 |. 



In the approximation when the EM coupling is taken to be constant in the cor- 
rection, the contribution of leptons is given by 

Kf (ml) = 3 ( 1 + l^) ln ^l = 3 (! + ^) • 1-98738 = 5.9725 . (52) 

The hadronic part is more involved. In the energy range from M T to m^ the number 
of active quarks is 4 or 9^ = 0. The partonic part of quark contribution reads 

. mi/ 9n 9 / ,3a\ mi /10 17 a \ ml , , 

Aff(™S) = W.^(l + *4-)b.4= ( T + Ys-Jm^l = 6-62896 (53) 



29 



where we use a = a. The result is independent of strong coupling constant (the parton 
model result without real QCD interaction). The quark part beyond the partonic 
approximation requires integration of the evolution trajectory of the strong coupling in 
nj = 4 effective theory. The initial value of the strong coupling is (M%) = 0.102001 
as was obtained from matching at M 2 for the c-quark contribution. In NLO one finds 
still a sizable contribution 

ASK) = (^-^-)4; ) = 0.54878. (54) 

3 54 7T 

The NNLO contribution proportional to a 2 s in eq. ([45|) 

ASK) = y^4?= 0-091848 (55) 



and the NNNLO contribution proportional to in eq. fl45|) 

A (% m 2) = ^200675 _ 335 } \ ( 3) = 63085J (3) = 0042699 (56) 
Tb K bJ V 23328 81 SV 7 Tb cb 

give only small corrections. Here 

A total correction to the parton model result (i.e. QCD contribution) 

AS^K) = A«K) + ASK) + A^K 2 ) = 0.644899 (57) 

is much smaller than the leading partonic result A^(m 2 ). For the EM coupling at 
m& one finds 

= 5^ "(5.9725 + 6.62896 + 0.644899) 

= 5<W- 13 ' 2464 ' (58 > 

Lepton and parton contributions dominate. Collecting all together one finds 

3n 3vT 



aW(m 2 b ) a( 4 )(M r 2 ) 
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13.2464 



And finally 



»7T »7T 

(32.7889 + 13.2464) = 46.0353 . (59) 

a a 



' = 132.152 (60) 



or 



a {A) (m 2 b ) = 1.037a. (61) 

This number can be used for the T-resonance physics analysis. 

Because decoupling is not explicit in mass-independent renormalization schemes 
there is another EM coupling parameter related to the scale m^. Upon changing the 
number of active quarks to rif — 5 one obtains 

Jkr^h +n2=ml0) - (62) 

The polarization function U b ^ ul1 (fj, 2 = ?7i 2 ,0) which gives the corresponding shift for 
the EM constant is written in terms of the effective strong coupling constant af\ml). 
A numerical value for af^{m 2 ) is obtained through matching the strong coupling at 
the scale mj. The running of the coupling a^\M 2 ) = 0.102 to m& = 4.8 GeV gives 
af\ml) = 0.06851 (a^(m 2 ) = 0.2152). Then matching at m b results in af\ml) = 
0.06869. With this number the result of matching for the EM constant is 

a = n 6/ui V = m 2 b , 0) = n Mr (/z 2 = m 2 , 0) + n Woop (/i 2 = m 2 , 0) 

= 0.00024 SM + 0.00358 ioop + 0.08587 + 0.03437 = 0.1241 . (63) 

The EM contribution is totally negligible. The loop contribution is rather small. The 
PT convergence of the direct contribution is not fast and is similar to the c-quark 
case. One has 

A fe (m 2 ) = 0.1241. (64) 
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Finally, the EM couplings of n/ = 4 and n/ = 5 effective theories in the vicinity of 
m 6 are related by 

^ AVD = -71^- 0.1241. (65) 



a®{ml) a^(m 2 b ) v 67 «( 4 )(mf) 
Explicitly one finds 

« (5) ^)= 132^138 = L000W4)(m ' ) - 

This difference can be safely neglected in applications for T-resonance physics. 

The uncertainty due to nib is tiny. Indeed, the error in the 6-quark mass leads to 
the uncertainty 

5A\ml) = -1(1 + af\ml)) 2 -^ = -\(l + afK))^ = ±0.030 . (66) 
3 nib 3 nib 

There are two reasons for such a smallness in comparison to the c-quark case: the 
electric charge of b quark \e b \ is two times smaller than |e c | and the relative uncertainty 
of the 6-quark mass mj (Snib/nib) is much smaller than that of the c-quark mass. 
Note that because the contribution A b (m 2 ) itself is small the relative uncertainty 
5 A b (ml) I A b (m 2 ) is huge. However, one cannot use it here. Even for A b (ml) = the 
uncertainty 5A b (m b l ) is basically 0.030. 



5.3 Running from to Mz- 



In this subsection we describe the evolution of the EM coupling constant cr 5 ) (ml) from 



mj = 4.8 GeV to Mz = 91.187 GeV. Various contributions according to eqs. ([45| fi8D 
are: 

a) Leptonic contribution 

3«\, Mr 



A[f = 3 ( 1 + ~) In = 17.6966 ; (67) 
\ An J ml 
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b) Leading quark partonic, a s -independent, contribution 



c) the NLO contribution with a^ 5 -*(m^) = 0.068694 as initial value for the evolution 



trajectory 



d) the NNLO contribution proportional to a 



3 144 

e) the NNNLO contribution proportional to a 3 



S = ~^= 0-UBlSi (70) 



s 



( 3) = Z257543 _ 620 x (3 ) = 6808gJ (3) = Q 2g 4 (n) 
bZ V 46656 81 sv 7 6Z bZ 

Here 

Jfi> = /?(«.?>(.)>"-• 

A total QCD correction to the partonic result 

A (hadcor) = A (l) + A (2) + ^(3) = ± ^ ^ 

is small compared to the leading quark partonic, a s -independent, contribution given 
in eq. (0). The total effect of running on the interval from mj to Mz 

A l b f + A[§ + A { b h z adcor) = 17.697 + 21.6048 + 1.1772 = 40.479 (73) 

is dominated by leptons and by the quark partonic contribution. We find the EM 
coupling at Mz expressed through the EM coupling at in the form 

= J? Tx ~ 40.479 . (74) 

This equation gives the relation between the running EM couplings necessary for 
applications in b- and Z-physics. 
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Collecting together eqs. P0|J58|J6^J7^ ) we find the absolute value of the running 
EM coupling at M z expressed through the fine structure constant a 

37T 37T 



2- = 32.7889(Wc/i M T ) - 13.2464(run m c 2m b 



and 



ct( 5 )(Mf) a 

3n 

- 0.12Al(match m b ) - 40.479(run m b 2M z ) = 86.6384 (75) 

a 



1 1 

-- 86.6384/(3tt) = 137.036 -9.1926 = 127.843. (76) 



a( 5 )(Mf) a 

The result can be written as a relation between the running EM coupling and the 
fine structure constant 

«( 5 )(M|) = 1.0719a. (77) 
This number can be used for the Z-boson peak analysis. 

6 Summary of results 

In this section we give a brief summary of the calculation paying attention to uncer- 
tainties of the results. 

The uncertainty for the low energy normalization value at the r mass is given 



in eqs. (|42| . |43|j4^ ) . It is largely dominated by the uncertainty due to the c-quark 



contribution. Adding uncertainties 

5A( 4 )(M T 2 ) = ±0.078^ ± 0.330, 

in quadrature one finds 

5A (4) (M r 2 ) = ±0.339 

and 

0.339/(3tt) = 0.036. 
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Range 


(0 - M T ) 


M T — m b 


m b -M z 


total 




22.011 


5.973 


17.697 


45.681 



Table 1: Leptonic contributions. 



Matching 


A uds (M T ) 


A c (M r ) 


A b (m b ) 


Value 


9.9312 ±0.078 


0.8468 ±0.330 


0.1241 ±0.030 



Table 2: Matching at different scales. 
Finally one obtains for the low energy normalization value the result of the form 

' 133.557 ±0.036. (78) 



ct( 4 )(M T 2 ) 

The total error is dominated by the uncertainty due to the c-quark matching contri- 
bution which is mainly given by the uncertainty of the c-quark mass. 

For the scales m& and Mz, the errors due to running, which are basically because 
of the uncertainty of the coupling constant a s , should be included. The running itself 
is precise because /3-functions in eqs. (^5| , |46|) are computed up to high order of PT and 
the coupling constant a s is rather well known. The dominant contribution comes from 
leptons and partonic quarks and is independent of the genuine QCD interaction (see 

2 

Tables The EM terms give a tiny correction. In Table f|the quantity J^l a s (s) ^ 
with running for a s (s) in different orders and with or without EM contribution to the 
strong /5-function is presented. The inclusion of EM terms slows down the decrease of 
the strong coupling and the integrals are slightly larger; still it is completely negligible 
numerically. In the leading order we have the uncertainty in integrals due to errors 
of the initial value of the strong coupling 

^ = TTTFI Sa - (79) 

i + a s p ij ab 
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Power 


M T — m& 


m b - M z 


total 


A(°) ~ a° s 


6.6290 


21.605 


28.234 


A (1) ~ a l s 


0.5488 


1.078 


1.627 


A( 2 ) ~ a* 


0.0918 


0.102 


0.194 


A (3) ~ a 3 s 


0.0043 


-0.003 


0.001 


total sum 


7.2739 


22.782 


30.056 


/\ (hadcor) 


0.6449 


1.177 


1.822 



Table 3: Running of powers of a s . 

with L a f, = ln^/yU 2 . This equation suffices for error estimates of the QCD contri- 
bution into running. NNLO and NNNLO give only small corrections. One can find 
uncertainty of the running by varying the initial value of a s . The numerical results 
are close to the estimate given in eq. (f79|). Eq. (|79|) , however, has an advantage of 
being analytical and simple that makes the error evaluation more transparent. 

At mb (and Mz) errors due to running and due to matching the light quark con- 
tribution at M T are not independent: both are determined mainly by the uncertainty 
in a s . Therefore these errors should be added linearly. 

For the interval from M T to mj one has from eq. (ff9|) 

SA hadcor^ = 4565as = Q.025 

and the total error (with linearly added errors for matching the light quark contribu- 
tion at M T and running) is 

6A<®(ml) = ±(0.078 + 0.025),^+™ ± 0.330 c ± 0.030 b 

= ±0.103 MflAt+Mm ± 0.330 c ± 0.030 fe . 
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Order 



LO 
NLO 
NNLO 
NNNLO 



J$ with EM 



0.169104 
0.165542 
0.164938 
0.164671 



1$ without EM 



0.169100 
0.165538 
0.164934 
0.164667 



Table 4: The quantity I^J from eq. (51) in different orders of strong /3-function and 
with or without the EM contribution to the strong /3-function. 

Adding independent errors in quadrature one has 

5A (5 \m 2 b ) = ±0.347 



and 



0.347/(3vr) = 0.0368 



Finally, one finds the uncertainty for the EM coupling at m b 



= 137.036 - (4.89766 ± 0.0368) = 132.138 ± 0.0368 

For the scale Mz the error due to running is estimated by 

11 



SA 



hadcor | 
tZ \a s 



8I^(a s (m b )) = 13.Q5a s = 0.074 



which leads to 



5A (5) (M|) = ±0.078^ ± 0.330 c ± 0.030 fc ± 0.074 r 



±0.152/,; 



light+run 



±0.330 c ± 0.030 6 . 



(81) 
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Adding independent errors in quadrature one has 

5A (5) (M|) = ±0.3646 

and 

0.3646/(3tt) = 0.03869 
These estimates give the error for the coupling at Mz 

127.843 ±0.039. (82) 

Eq. (|52]) is a final result. However, it cannot be directly compared with the results of 
the standard analyses because the quantity in eq. (jS2]) is defined in a different scheme. 
We consider the uncertainty in eq. ([82]) (the part of which related to the running is 
estimated analytically for transparency) as rather conservative. 



7 Comparison with other schemes 



With the number from eq. (|8~2"D one can find the on-shell parameter a os at Mz used 
in the literature. Indeed, because of the relation 

- + II os (q 2 ) = + n( 5 V, q 2 ) . (83) 



a os (q 2 ) a "'" a^ 5 \fi 

one has to compute n^ 5 ^(/x 2 ,q 2 ) in rif = 5 effective theory at the point q 2 ~ Mf. 
Note that we have restored a notation q 2 in addition to the usual quantity q 2 : the 
new variable q 2 will be used in Minkowskian domain. For computing the leading part 
of U^(fi 2 , q 2 ) in the kinematical range fj 2 ~ q 2 ~ M| one can consider all five active 
quarks (u, d, s, c, b) and all three leptons as massless and use eq. (|T0|) with the only 
change because of a different number of active quarks which is now 5 instead of 3. 
This change affects only 0(a 2 ) order in eq. (|I~0"D and changes nothing for leptons in 
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the NLO approximation. One has for the generic light polarization function of quarks 
n«'*-V, q 2 ) = In ft + | + « s L £ + g - 4<(3)) 

+al W-t- 2 + (_ - -n, - 4A( B/ )C(S) ) In jL 

41927 _ 3701 _ ,829 _ 19 N 25 | 

864 1296 ' V 18 9 V W 3 SV V V ; 

with Po(nj) = (11 — 2/3n/)/4. For a more accurate evaluation of IT 5 )(// 2 ,q 2 ) at the 
scale we retain the leading corrections due to the b-, c-quark and r-lepton masses 
and the leading correction due to the top quark contribution. One finds 

II (5) (/i 2 ,q 2 ) = m li9ht - lept (ii 2 ,q 2 ) + —u li9ht - quark (fi 2 , q 2 ) 

3 

16m? 4 6m 2 6M 2 . .-s. 9 9 . 

- o^r - o^f r + Awn (5) (/i 2 , q 2 ) • (85) 

3 q^ 3 q^ q z 

Note that the power correction due to a quark mass is exact up to 0(a 2 ) order when 
expressed through the pole mass. The corrections due to the top quark contribution 
for the quantity Il^ 5 ^(/i 2 , q 2 ) at q 2 » Mf can be computed as a power series in q 2 /m|; 
the expansion parameter q 2 /m| is small at the point q 2 = Mf for m t = 175 GeV. 
Indeed, retaining only the leading term and first corrections one has 

A w Il( 5 )(M 2 ,q 2 ) = -±4 I 1 + ^ ] (M 2 Z ) - ^4} ■ (86) 

15 mf [ 81 28 mf J 

A typical expansion in eq. fl5S|) reads 

A (t) n (5) (M|, M|) = -0.0724 - 0.0138 as + 0.0021 Mz = -0.0841 (87) 

with obvious notation indicating the origin of different contributions. We do not take 
into account bosons therefore the VT-boson loops should be analyzed separately. 

To calculate the on-shell coupling a os (q 2 ) at the scale Mz using eq. (|8~3"D one can 
use either q 2 = Mf (Euclidean definition) |50| or q 2 = — q 2 = Mf with taking the real 
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part of the correlator (Minkowskian definition). The Minkowskian definition is usually 
discussed in the literature. Note that we calculate not the e + e~ scattering amplitude 
at the total energy Mz (q 2 = Mf ) which definitely should be taken at a physical 
point on the cut in the case of cross section calculations, but the coupling constant 
which parameterizes this amplitude at the scale Mz- Within the RG approach the 
scale of the parameters of an effective theory valid in a given energy range should not 
coincide with any actual physical value of the energy or momentum squared (see e.g. 



First we use a Euclidean definition for the on-shell coupling which is consistently 
perturbative and requires computation of the correlator IT( 5 )(Mf,q 2 ) in a deep Eu- 
clidean domain for q 2 = M§. Using eqs. ( |35"| , |8~6"D one finds an expansion 



II (5) (Mf , Mf ) = 11.1111 - 0.03097 as + 0.00112 a 2 

-0.00168.ea/ - 0.00554 - 0.00304 c - 0.00228 T - 0.0841 t 

= 11.0796 - 0.01086 6cr - 0.0841 4 . (88) 
Note that the EM correction is numerically of the order a 2 . Still these corrections are 



very small. From eq. (j83|) we find 

37T 37T 



+ n {5) (M§,Mf) 



a os (M 2 z ) «( 5 )(Mf) 

= — - 86.6384 + 11.0796 - 0.01086 6cr - 0.0841 t 

a 

= -75.5588 - 0.01086 bcr - 0.0841 4 . 

For clarity we retain the contribution of power corrections separately for further 
comparison with the results in Minkowskian domain. One has numerically 

1 - - (75.5588 + 0.01086 fecr + 0.0841 t )/(3vr) 



a os {M 2 z ) a 

= 137.036 - 8.0271 = 129.009 . 
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Because the error estimate in eq. ([82]) is not affected by this change of scheme (too 
small and rather precise contributions are added), the final result for the on-shell 
coupling within Euclidean definition reads 

y-r-ppr = 129.009 ± 0.039 . (89) 

However, the reference values for the on-shell coupling available in the literature are 
given in Minkowskian domain for q 2 = — q 2 = M§, i.e. for the real part of the 
correlator U^(M^,q 2 ) computed on the physical cut. Within the approximation 
used, going to the Minkowskian domain of momenta q 2 changes only the a 2 order 
term and power corrections in eq. (|8"8|) . Indeed, in eq. ( |S~4"D the only term which is 
numerically affected by the transition to the Minkowskian domain is ln 2 (yU 2 /q 2 ) with 
the following change as compared to the Euclidean result 

Instead of eq. ([88]) one finds 

Re {n (5) (M|,-M|)} = 11.0796 
- 0.04893^2 + 0.01086 6crM mfc + 0.08828 tMmfc = 11.1298 (91) 



and 



Re z — ■ — ^ = 137.036 - (86.6384 - 11.1298)/(3tt) = 129.024 . 

a os (-M|) 



The final result for the on-shell coupling with Minkowskian definition is 

Re \^ = 129.024 ± 0.039 . (92) 

«os(-M 2 ) v ; 

The difference between the central values for the couplings in eq. ( |8"9"P which gives the 
Euclidean definition and in eq. (|92]) which gives the Minkowskian definition is 0.015. 



Note that the Euclidean definition was considered in ref. 50 where the numerical 
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difference about 0.02 from the Minkowskian definition has been found from rather a 
simplified treatment. It is close to the present more accurate result 0.015. Note that 
the point Mz is safe for the PT calculation in Minkowskian domain for the approx- 
imation used (no singularities of the spectrum near this point). At other points it 
is not so even in the approximation we work. For instance, if the T-resonance mass 
mx is taken as a reference scale then the Euclidean definition is equally applicable at 
this point being still perturbative while the Minkowskian definition faces the prob- 
lems that the polarization function on the cut is not smooth. A phenomenological 
approach based on direct integration of data fails because of the fast change of the 
spectrum at the location of the T resonance which makes the integration with prin- 
cipal value prescription for regularizing the singularities ill-defined. A theory based 
approach within PT fails at the point mr because PT calculations for the correlator 
near the threshold on the physical cut (mx ~ 2m&) are not reliable. Therefore, the 
Minkowskian definition is not uniformly applicable at every scale. 

The present paper result given in eq. ( [92]) differs from some recent determinations 
based on the use of experimental data for integration over the low energy region. For 
the result of ref. [|16j 

Re y- — ^ = 128.925 ± 0.056 (93) 

the number obtained in the present paper and given in eq. (|92| ) almost touches the 
reference value eq. (|93| ) within la (a is a standard deviation). The results of some 



other groups are concentrated around the same central value as in eq. (93) but with 



essentially smaller errors. For further comparison, we use the result of ref. jl7| 



Re — ^ = 128.927 ± 0.023 . (94) 

a os (-M 2 z ) K ! 



The difference between the value from eq. (^) and the reference result in eq. fl9~4] ) is 
129.024 — 128.927 = 0.097 which constitutes 2-4<r and can be significant. Therefore 
we discuss the difference in more detail. 
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The usual parameterization of the fermionic contributions to the on-shell running 
EM coupling at M z reads 

Re t^—tttt = — ( 1 - Aai op - Aa£d - Aa top 



a os {-M 2 z ) a 

(note that 



Re U*) * ' 



y a os (-Ml)J ^ Rea os (-M z ) 
though the difference is tiny). The total leptonic contribution to the electromagnetic 
coupling constant at M z given in the last column of Table [1] reads 

A lep (M 2 z ) = 45.681. 

The leptonic part of lT 5 )(M§, -Mf) reads 

Re {u^ lep {M 2 z , -Mf)} = 3Re (n^'-^(Mf , -M|)) + 



4.9988 + 0.0023 = 5.0011. 



The leading order contribution is equal to 5 while the EM and r-lepton mass cor- 
rections are small. For the total leptonic contribution to the on-shell coupling in the 
Minkowskian domain one finds 

Aaiep = ^(45.681 - 5.001) = 314.974 x 10~ 4 
on 



which is close to the number of ref. [17] 



Aaiep| r ef = (314.19 + 0.78) x 10" 4 = 314.97 x 10" 4 . 

For the top contribution in the Minkowskian domain one finds from eq. ( jS6|) (see also 
eq. d?D) 

Aa top = ^(-0.0883) = -0.68 x 10" 4 
on 
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while the number of ref. [17) with a more accurate account for the higher order 
corrections is 

Aa top | ref = -0.70x nr 4 . 



The difference is small and is neglected. From the numerical value given in eq. ( |9~2"D 
the total contribution of fermions into the shift of the EM coupling is determined to 
be 

Aai cp + Aa top + Aafi = 1 - a ■ (129.024 ± 0.039) = 0.0584664 ± 0.000285 . 

Taking Aai ep and A« top as exact quantities (no errors) one finds the following nu- 
merical value for Aa h ^ d 

Acfi = (0.0584664 ±0.000285) - 0.031497 + 0.000068 

= (270.37 ± 2.85) x KT 4 (95) 



while the result of ref. |17J is 

„(5) 



Aa^|i = (277.45 ± 1.68) x KT 4 , (96) 



and the number of ref. [[HJ is 

Aa£dh = (277.6 ± 4.1) x 10" 4 . (97) 



The difference between the central values given in eq. fl95|) and eqs. (|96| , |97|) is about 
2-4a depending on the numerical value of the error quoted 

277.45 - 270.37 = 7.08 » 2.5 • 2.85 « 4.3 ■ 1.68 » 1.7 ■ 4.1 . 

This difference can be significant. Therefore, we discuss the sensitivity of our predic- 
tion fl9"5] ) (and of eq. ( |S2"D from which it is uniquely obtained) to the numerical values 
of parameters used in the theoretical calculation of the present paper. If m s = then 
the p- and ^-channels should be theoretically degenerate because there is no reason 
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for the difference. This means that besides vanishing explicit corrections due to m? s 
in eq. ( p2[) one should identify m^ s with the resonance in the nonstrange channel, i.e. 
numerically substitute = m p into the solution for the IR modifying parameters 
in eqs. (|20|) . With such changes one finds the result for A uds (M 2 ) in the form 

A uds (M 2 T )\ ms=0 = 10.23 (98) 

that generates a numerical shift about 0.3 in the value of A uds (M 2 ) as compared to 
the result for nonvanishing strange quark mass in eq. (|30|). Note that if the direct 
integration of low energy data is used then the full dependence of the results on m s 
is lost. Only the PT high energy tail depends explicitly on m s but this dependence 
is weak. The change in Aa[ 5 ^ corresponding to the result in eq. (|98|) is 



A«£dU fl =o - Aa® = 0.3^- = 0.000232 = 2.3 x 10" 4 . 

on 

The use of the numerical value m c = 1.6 GeV for the c-quark mass instead of m c = 
1.777 GeV generates the 0.33 shift in the value of quantity A C (M^) that leads to the 
following change of Aa h ^ d 

AcfiL^i.e - Aafi = 0.33^- = 0.000256 = 2.6 x 10" 4 . 

Note that this change cannot be found if the direct integration of contributions of 
actual charmonium resonances is performed. The total shift due to such a change in 
these parameters compared to eq. (|9~5|) is 

A4adU=o,m c =i.6 = (270.37 + 2.3 + 2.6) x 10~ 4 

= (275.27 + 2.85) x 10~ 4 . (99) 

The change of the numerical value of the strong coupling constant a s from a s (M 2 ) = 
0.318 to a s (M T 2 ) = 0.335 gives a 0.152 shift in A (5) (Mf) (according to our error 
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estimates in eq. (|8lD) to end up with 

Aa^i| m , = o,m e =i.6,a.=o.335 = (270.37 + 2.3 + 2.6 + 1.2) x 1(T 4 

= (276.47 ± 2.85) x 10" 4 . (100) 

This result agrees with other estimates within la. The set of numerical values for the 
relevant parameters used in eq. (|100|) is rather close to the set used for obtaining the 



value in eq. (g5J) (m s = 130 MeV, m c = 1.777 GeV, af(M r 2 ) = 0.318). The total shift 
in Aa^ad f° r the new se ^ °f parameters in eq. (|100|) is larger than the total error given 



in eq. (95) because the total error is computed in quadrature and the change of the 
spectrum due to m s = (which makes all three light quark channels degenerate) has 
not been included into the total error. To definitely distinguish between the results 



of eq. (|100|) and eq. (|95|) more precise numerical values of parameters are necessary. 



Within the present paper approach we use virtually no real data on cross sections 
but rely on the numerical values of several theoretical parameters which are extracted 
from such data. These parameters are the strong coupling constant, quark masses, 
vacuum condensates. It is generally believed that the real data can be properly 
described with these parameters if theoretical formulas are sufficiently accurate. In 
the case of computing the hadronic contribution into the photon vacuum polarization 
function in Euclidean domain a theoretical description is pretty accurate because PT 
is applicable and very precise - in fact, the PT results in this area are almost the 
best ones available among all PT calculations. An additional reason for such a high 
theoretical precision is that for the calculation of Aa^ ad in Euclidean domain one 
extracts only very general information encoded in the data - just the integral over 
the entire spectrum with a smooth weight function and no details of the behavior 
over specific energy regions. This is the situation when global duality, which is exact 
by definition (hadron and quark descriptions are supposed to be exactly equivalent in 
principle), is applicable and is under a strict control numerically within PT. However, 
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our calculation shows that at the present level of precision the result for Act had is 
rather sensitive to the numerical values of the parameters m c and a s which should 
be fixed from the data. The uncertainties of these parameters can be reduced both 
with better data and better theoretical formulas for extracting the numerical values for 
these parameters from the data (especially m c ) while the theoretical framework for the 
calculation of the hadronic contribution itself is already very precise. Using the result 
given in eq. fl9~5|) and formulas for radiative corrections to the Weinberg angle from 
ref. |S2j (assuming Minkowskian definition for the on-shell coupling) we find that the 
central value of the Higgs-boson mass moves from the reference value Mh = 100 GeV 
for Aafi = 280.0 x 10" 4 to M H = 191 GeV for the value Aa^ = 270.37 x 1(T 4 
found in the present paper. 

8 Conclusion 

The technique of calculating Aa h ^ d within dimensional regularization and minimal 
subtraction is straightforward in PT. It heavily uses the renormalization group which 
is the most powerful tool of modern high precision analyses in particle phenomenology. 
Because PT is not applicable only at low energies one should modify only the IR region 
of integration for light quarks: a numerical integration of data at energies higher than 
2-^3 GeV is equivalent to the theoretical calculation in PT if both data and theory are 
properly treated. The present calculation uses virtually no explicit scattering data but 
the values of the lowest resonance masses for the light-quark vector channels. Other 
experimental information is encoded through the numerical values of the coupling 
constant, quark masses, some vacuum condensates. 

Minkowskian definition of the on-shell coupling constant is deficient and not ap- 
plicable at some points. Both MS and on-shell coupling constants in the Euclidean 
domain can be reliable determined with the use of theoretical formulas already es- 
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tablished in high orders of PT. In view of future high precision tests of SM at Mz 
and two-loop calculations for the observables in this region it seems that the param- 
eterization of the theory with the running EM coupling in the MS scheme is most 
promising. 

The main uncertainty of the hadronic contribution into the running EM coupling 
constant at Mz comes from the error of the numerical value of the c-quark mass 
m c . The uncertainty of a s is less important. Unfortunately, the c-quark mass is a 
quantity which is very complicated to study. The reason for that is its numerical 
value close to the strong interaction scale of the order of p-meson or proton mass. 
Therefore, m c should be treated exactly in theoretical formulas, almost no simplifying 
approximation is applicable in the kinematical range of energies of order m c . The 
presence of a mass usually makes the loop calculations within PT technically more 
difficult. Near the cc production threshold where the mass is essential and where 
its numerical value can be reliably extracted from accurate experimental data, the 
Coulomb interaction is enhanced that requires to take it into account exactly while the 
c-quark mass is too small for NRQCD to work well. Finally, the nonPT corrections 
due to vacuum condensates within OPE are essential numerically in this energy range 
but they are not well known because they are given by the gluonic operators ||53| . 
And though the coefficient functions of the relevant operators up to dimension eight 



are calculated ||54|| , the numerical values of their vacuum condensates are poorly 
known. These reasons make the accurate determination of the c-quark mass difficult. 
The uncertainty related to the contribution of the c quark to the hadronic vacuum 
polarization is additionally enhanced because of the c-quark large electric charge. 
For the b quark, for instance, all the above problems are much less severe. Therefore, 
the c-quark physics plays an essential part in the Higgs search through radiative 
corrections. 
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